Multisetting Bell-type inequalities for detecting genuine tripartite entanglement 
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In a recent paper, Bancal et al. put forward the concept of device-independent witnesses of 
genuine multipartite entanglement. These witnesses are capable of verifying genuine multipartite 
entanglement produced in a lab without resorting to any knowledge of the dimension of the state 
space or of the specific form of the measurement operators. As a by-product they found a three- 
party three-setting Bell inequality which enables to detect genuine tripartite entanglement in a 
noisy 3-qubit Greenberger-Horne-Zeilinger (GHZ) state for visibilities as low as 2/3 in a device- 
independent way. In this paper, we generalize this inequality to an arbitrary number of settings, 
demonstrating a threshold visibility of 2/7T ~ 0.6366 for number of settings going to infinity. We 
also present a pseudo-telepathy Bell inequality achieving the same threshold value. We argue that 
our device-independent witnesses are optimal in the sense that the above value cannot be beaten 
with three-party-correlation Bell inequalities. 



I. INTRODUCTION 

Quantum theory allows correlations between remote 
systems, which are fundamentally different from classical 
correlations Quantum entanglement is in the heart 
of this phenomenon Q. Already two entangled particles 
give rise to correlations not reproducible within any local 
realistic theory Q. However, moving to more particles 
a much richer structure and various types of entangle- 
ment arise [J] suggesting novel applications such as quan- 
tum computation using cluster states @, sub-shotnoise 
metrology @, or multiparty quantum networking Q- In 
these tasks, genuinely entangled particles offer enhanced 
performance. Hence, it is a central problem to decide 
whether in an actual experiment genuinely multipartite 
entanglement has been produced, or alternatively, the 
entangled state prepared in the laboratory could be ex- 
plained without requiring the interaction of all particles. 
In the latter case, we say that the state created is bisep- 
arable. Focussing on the tripartite case, a biseparable 
state pbs can be written as 



Pbs 



^2,Pi\<t>i)(<t>i 



(1) 



where the pure states <pi are separable with respect to one 
of the three bipartitions 1 1 23 , 12 1 3, 13| 2, and the weights 
Pi > add up to 1. For more than three parties, the 
generalization is straightforward. 

Several experiments have been conducted so far gen- 
crating multipartite entangled photonic states up to six 
photons (for instance, Ref. 3 generated a Dicke state 
of six photons). One of the traditional approaches to 
decide on the existence of genuine multipartite entan- 
glement consists in performing a complete state tomog- 
raphy, and then deducing the kind of entanglement di- 
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rectly from the density matrix using witness operators. 
Alternatively, the experimentalist may measure cleverly 
chosen witness operators, thereby reducing the number 
of correlation terms to be measured in the actual ex- 
periment Q. However, a common drawback is that in 
both cases the experimentalist needs to have a precise 
control over the system on which the measurements are 
performed. 

Remarkably, there is another route, avoiding the above 
problem, building on the seminal work of John Bell [lj: 
Bell expressions are linear functions of joint correlations 
enabling one to say important things in a black box sce- 
nario about the dimension of the systems, the states in- 
volved, or the kind of measurements performed. In par- 
ticular, it is possible to decide on the presence of gen- 
uine multipartite entanglement based on merely statisti- 
cal data (that is, without relying on any knowledge of the 
implementation of the devices involved in the measure- 
ment process) [l(J [Hj]: if a Bell value, coming from the 
statistics of a Bell experiment, is bigger than a certain 
value achievable with measurements acting on bisepara- 
ble quantum states, then we can be sure that the state 
in question is genuinely multipartite entangled. This ap- 
proach has been formalized more recently by Bancal et 
al. [111 ], coining the term device- independent witnesses 
of genuine multipartite entanglement witnesses for such 
Bell expressions (for more details we refer the reader to 
that paper). 

As a simplest illustration of a device-independent wit- 
ness of genuine tripartite entanglement, let us represent 
the Mermin polynomial 12] in terms of three-party cor- 
relators, 

h =(A ® B ® C ) ~ (in Ci) - (Ai <g> -B ® Ci) 

-(Ai®Bi®C ), (2) 

where (A ce <g)-Ba(g>C' 7 ) designate the expected value of the 
product of three ±l-observables, A a , Bp, C 1 . It has been 
shown in Ref. [l(|, that I2 < 2\/2 for biseparable quan- 
tum states (B2 = 2v2)j whereas the maximum quantum 
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value saturates the algebraic limit of 4 (Q2 = 4), hence 
the violation of the bound B2 implies genuine tripartite 
entanglement. Note that this reasoning holds true inde- 
pendently on the size of the Hilbert space dimension or 
on the type of measurements carried out. Hence, Mer- 
min inequality serves as a device- independent witness of 
genuine tripartite entanglement [llj . Let us now take the 
noisy 3-qubit GHZ state, 

p(V) = V\GHZ){GHZ\ + (l-V)± (3) 

where \GHZ) = (|000) + jlll>) /y/2 is the 3-qubit GHZ 
state [13J, and V is the visibility parameter. The mea- 
surements achieving the bounds Q2 and B2 correspond 
to traceless observables, entailing the threshold visibility 

V = B2/Q2 = l/v2< Hence, genuine tripartite entan- 
glement in the noisy GHZ state (for V > l/y/2) can be 
detected in a device-independent way. 

More recently, however, Bancal et al. [ll| managed to 
lower the threshold visibility of the noisy 3-party GHZ 
state to V = 2/3 by considering a three-party three- 
setting Bell inequality, which can be considered as a 
three-setting generalization of the two-setting Mermin in- 
equality. Note that similarly to the Mermin inequality, 
the Bancal et al. inequality extends to more than three 
parties as well 

In the present paper, we generalize the three-setting 
three-party Bancal et al. inequality to an arbitrary num- 
ber of settings m, exhibiting the threshold visibility 

V = l/(msin(7r/2m)), which approaches V = 2/ir for 
large number of settings. This generalization is discussed 
in section [Hi whereas another family of Bell inequalities, 
based on the extended parity game [14|, is discussed in 
section Hill Notably, this game exhibits pseudo-telepathy 
[l5j , and the corresponding Bell inequality has the same 
performance (for m a power of 2) as our Bell inequality 
of section |TTJ 



We next exhibit a lower bound on the quantum max- 
imum, Q l m — m 3 /2, as a function of number of set- 
tings m. Then an upper bound is given on the bisep- 
arable quantum maximum, which is shown to be at- 
tained by von Neumann-type projective measurements, 
B m = m 2 /(2 sin(7r/2m)). This implies the threshold vis- 
ibility V = B m /Q m < B m /Q l m = l/(msin(7r/2m)) tend- 
ing to 2 /tt in the limit of large number of measurement 
settings. 

We wish to note that Bancal et al. (Appendix C in [Til ]) 
presented a biseparable model simulating all the single- 
party expectations (A), (B), (C), two-party correlators 
(A <E> B), (A® C), (B <E> C) and three-party correlators 
(A(g)B(g>C) , achievable with von Neumann measurements 
on the noisy 3-qubit GHZ state Q of visibility V < 1/2. 
Within this biseparable model all of the three parties may 
share local random variables, but at most two parties can 
share a quantum state at a given time. It can be shown 
that if we are content with simulating only the three- 
party correlators, then the threshold visibility becomes 
a higher value, V = 2/ir. This implies that it is not 
possible to detect genuine tripartite entanglement in the 
3-qubit GHZ state in the range V < 2/ir by applying von 
Neumann measurements and considering Bell expressions 
which are sums of three-party correlators. In this sense, 
our family of Bell inequalities is optimal, giving V — > 2/n 
when m goes to infinity. 



Lower bound on the quantum maximum, Q} m . If 
each of the participants performs a von Neumann pro- 
jective measurement on one component of a shared 3- 
qubit GHZ state, the tripartite correlation of their mea- 
surement values can be written as (see for instance Ap- 
pendix C in Ref. [ll|): 



II. MULTISETTING TRIPARTITE 
BELL-TYPE INEQUALITIES 

Let us introduce the m-setting tripartite Bell expres- 
sion, 



{A a ^Bp^C 1 ) = sia0£ sin 9f sin 0% cos{<p£ + tp% + <p°), 

(6) 

where the ath, /3th and 7th measurement operator A a , 
Bp and C 1 of Alice, Bob and Cecil, respectively are given 
as: 



I m = ^2 M al3j(A a ® Bp 
a,/3,7=0 



a 



7/' 



where the matrix of Bell coefficients is defined by 



M a p 1 — cos 



■(a + /3 + j- A) 



(4) 



(5) 



where indices a, and 7 may take the values of 
0,1,..., m — 1, and A may be any real number. By 
choosing m = 2 and A = 0, the Mermin polynomial ([2]) is 
recovered. On the other hand, for m = 3 and A = —1/2, 
we obtain the polynomial of Bancal et al. [Tl| (apart from 
an irrelevant multiplicative factor). 



A a = cos sin 0^ a x + sin ip£ sin 0^a y + cos 9^a z 
Bp = cos <fp sin 9p a x + sin ipp 3 sin dp a y + cos 9p a z 
C 1 = cos tp^j sin 9^a x + sin tp^ sin 9^a y + cos 9^a z , (7) 



where a x , a y and a z are the Pauli operators. 



With the choice of Q&=Q? = VP- 

M H- h L 



and = <p» = 

ip"^ = 7r(/z — A/3)/m each tripartite correlation will take 
the same value as the Bell coefficient to be multiplied 
with, and the quantum value of the Bell expression will 
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be easy to calculate: 



m — 1 m — 1 m— 1 



M. 



{A} 



a=0 /3=0 7=0 
m— 1 m— 1 'm— 1 

EEE cos 

a=0 £=0 7=0 
1 m— 1 m— 1 m — 1 ✓ 

jEEEf- 

a=0 0=0 7=0 

to 
~2~' 



(a + /3 + 7- A) 



2tt 

— (a + /3 + 7 - A) 

TO 



(8) 



This value Qj n is a lower bound for the maximum quan- 
tum value Q m . 

For the maximum of the biseparable value first we will 
give an upper bound (£>^J, then we will prove that this 
bound can be saturated, that is, (B l m = B^ = B m ). 

Upper bound on the biseparable quantum value, 

B^. The value to be calculated is, 

m — 1 m — 1 rn — 1 

B m = max ^ X X M a/3l A a (Bp ® C* 7 ) 

Q =0 0=0 7=0 
m — 1 m— 1 

= max ^ J] M if <A> ® ^r). ( 9 ) 

/3=0 7=0 

where we used the fact that Bell inequality Q is linear in 
the correlators, and that a biseparable density matrix (fT]) 
is a convex combination of pure states, hence it is enough 
to take the 3-party correlators in the form (A a ) (Bp^Cj) . 
In Eq. © each of A a may take the value of either +1 
or — 1, Bob and Cecil may share any quantum state and 
perform measurements on them, the operators of their 
measurement settings are Bp and C 7 , respectively, and 
the coefficients of the two-partite Bell inequality, which 
depends on the actual choice of A a are: 



'(/3+l)m ~ ~ Mpf^, that is the last element of each row 
is the same as minus one times the first element of the 
preceding row. These properties are v ery similar to the 
properties defining circulant matrices [17j j . whose eigen- 
vectors are independent of the actual values of its ele- 
ments, and therefore whose eigenvalues are very easy to 
derive. There are just two differences. In the case of 
the circulant matrices the elements are shifted not to the 
left, but to the right. Furthermore, they do it cyclically, 
that is there is no change of sign when the last element 
takes the first place in the next row. The first difference 
is easily corrected if we rearrange Cecil's measurement 
settings into the opposite order. Fortunately, the change 
of sign of the matrix element poses no serious problem 
either, because it can be shown that the eigenvectors of 
these modified circulant matrices are also independent of 
the actual values of the elements of the matrix, they are 
given as: 



(11) 



where j = 0, . . . , m— 1. The difference from the circulant 
case [17| is the 1/2 term in the exponent. To calculate 
the eigenvalues we only need the first row of the matrix. 
Therefore, we get the upper bound for the biseparable 
value as: 



B™ = max m 



m— 1 m—1 

53 5Z ^ a C0S 

q=0 7=0 



7r(a — A' — 7) 



u 1 



(12) 

where we introduced the notation A' = A — m + 1 . By 
substituting the ojj from Eq. pT]) and using identity 



M 



{A} 



5> Q 

«=0 



cos 



— (a + ^ + 7- A) 

TO 



(10) 



We note that due to the symmetry of the Bell expression 
under party exchange, it is enough to consider the case 
when it is Alice who may not share an entangled quantum 
object with the others. 

From the work of Ref. r JjJ it easily follows that for bi- 
partite correlation type Bell inequalities with an equal 
number of measurement settings per party, an upper 
bound for the maximum quantum value is the largest of 
the singular values of the matrix defined by the Bell co- 
efficients multiplied by the number of measurement set- 
tings. In the present case the matrix depends on the 



sum of its indices. Therefore, M, 



{A} 



M 



{A} 



that 



'(/3 + l) 7 ~ 1 "/3(7+l) 

is each row contains the elements of the preceding row, 
shifted to the left. From Eq. (fTU|) it is also clear, that 



nia — A — 7) 7r(a — A ) ir^/ 
cos = cos cos — 

TO TO TO 

7r(a — A') 7T7 
+ sm sin — 

TO 771 



we arrive at: 



B!t = max to 



a=0 



*(«- A '),ry_ L .-nJ 



TO 

ir(a - A' 



(13) 



rn 



(D{ + iDi) 
(Di+iDi) 



(14) 
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where 



Then we get 



7=0 
m—l 

4=E 



cos — COS 
to 



7=0 
m—1 



cos — sin 

to 



sin — cos 

TO 



^ = E 

7=0 
m— ] 

^ 4= y sin _i 



2n(j 


+ 1/2)7 




m 


Mi 


+ 1/2)7 




m 


Mi 


+ 1/2)7 




m 


Mj 


+ 1/2)7 



7=0 



in 



(15) 



However, 



n*+n> V- ■ 27r(j + 1/2 ±1/2)7 „ 

D 2± D 3= 2^ Sm Z = °' (I 6 ) 



7=0 



therefore Z?2 = Dl = 0, and 



TO 



m—1 

E 

7=0 



cos 2,( J + l/2 T l/2 h 



from which it follows that D\ = D 3 4 = for 1 < j < m— 2, 
D° = Dl= to/2, and Z?" 1-1 = -D™" 1 = to/2. To get 
the maximum we must take either j = or j — = m — 1. 
For j = we get: 



£ m = max — 



a=0 



A a e 



(18) 



If we have taken j = m — 1 instead of j = 0, we would 
have got the complex conjugate of the numbers whose 
absolute value has to be taken, which would have given 
the same result. In Eq. (fT8|) we have to add to vectors 
on the complex plane, each pointing towards corners of 
a regular polygon of 2to sides, and then we have to take 
the length of this vector. Each vector lies on a different 
diagonal of the polygon, but may point towards either 
direction depending on the value of A a . It can be shown 
that we get the largest value if the vectors taken in some 
order point towards consecutive corners. All such ar- 
rangements give obviously the same result. We get one 
of those arrangements if we take A a — 1. The result does 
not depend on A', as changing A' means only an overall 
rotation of the arrangement. Let us take A' = —1/2. 
Then the set of numbers will be symmetric with respect 
to the imaginary axis, therefore the real part of the sum 
will be zero, while the imaginary part will be positive. 



to 2 ^ . 7r(q + l/2) 

a=0 



2 sin 



9 m — i 

m ( 

■in ^- ^ V 



7r(a + l/2) . 7T 



sm ■ 



2m a= Q 

n(a + 1/2) 7T \ 
— ■ — - rns 



2m 



COS ■ 
to 2m/ 



9 m — 1 
TO x ^ 7Ta 
COS 



"2sin^ ^ 

2m a= Q 



2 sin J- ' 

2m 



(19) 



Here we have used that J2™=o C0S [ 7r ( Q! + 1/2)/to] = 0, 
and that cos(7ra/TO) = — cos[7t(to — a)/m\. 

Now we will show that this upper bound can be satu- 
rated. 

Lower bound on the biseparable quantum value, 



If S™=o S™=o MpyCp ■ C 7 is a certain number, 
where Cp and C 7 are Euclidean unit vectors, then there 
exist measurement operators giving the same number as 
the quantum value of the bipartite correlation type Bell 
inequality of coefficients M ^ 7 , applied on the maximally 
entangled state [lSj . In case of two dimensional vectors 



{^4} 



pairs of real qubits are sufficient. Let Mg 7 = 
all A a = +1 (see Eq. (fTU|) ). let Cecil's vectors be 

. . . , c _2 

Cry = 



with 



(20) 



and let us choose Bp optimally, that is Bp = 

J2™=o MpjCj/l J2™=o Mp^C^. Then the corresponding 
quantum value is: 



0=0 
m—1 

= E 

/3=0 



m—1 

E ^P-y C f 

7=0 

m—1 m—1 

E E cos 

a=0 7=0 



7r(a + P + 7 - A) / cos ^ 



in 



sm 



77".. 



(21) 



Now if we follow analogous steps to the ones we have 
taken calculating the value of B u , we will arrive at the 
same result. We can also easily see this if we compare 
Eq. (|2"Tj) to Eq. (jl"2"j) . The maximum value of the latter 
expression has been attained with A a = 1 and j = 0. By 
substituting these values, and also ojo from Eq. (jlip . we 
get almost the same formula as Eq. (|2"Tj) . indeed. The Uq 
complex numbers correspond to the same vectors on the 
complex plane as the two dimensional vectors appearing 
in Eq. (f2"Tj) . From the calculation of B u it turns out that 
the value does not depend on A', so the absolute value 
in Eq. (|2"Tj) does not depend on A — (3 either, therefore, 
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we may replace the summation in terms of (3 for a mul- 
tiplicative factor of m. The only remaining difference is 
the opposite sign of 7 in the cosine, but that will not 
affect the result either. 

As the upper and lower bound for the biseparable case 
are equal, the biseparable value itself is given by Eq. (|T9| . 
For the quantum value we have only proven a lower 
bound (see Eq. ©). Therefore, for this family of Bell 
inequalities the ratio of the quantum and the biseparable 
values (which equals the visibility threshold) satisfies: 



V = 



B„ 



< 



1 



(22) 



We believe that the lower bound Q l m we have given in © 
is actually the quantum maximum itself, and the above 
expression is valid as an equality. 



III. BELL-TYPE INEQUALITIES BASED ON 
THE EXTENDED PARITY GAME 



Now we define another family of multisetting tripartite 
inequalities giving the same ratio of B m j Q m as the right 
hand side of Eq. (|22|) . at least when m is a power of 2. 

The Bell coefficients may only take the values of zero, 
one and minus one, namely: 



if (a + + 7) mod ra/0, 
Mafij — { 1 if (a + /? + 7) /m is even, 
— 1 if (a + (3 + 7)/m is odd, 



(23) 



and a, /3, 7 = 1, . . . , m — 1. These Bell coefficients corre- 
spond to the so-called extended parity game considered 
in Ref. 14]. An equivalent definition, more similar to the 
definition of the Bell inequality ([5]) treated in section HTl is 
that M Q( 3 7 = cos[7r(a + /3 + 7)/m], whenever the absolute 
value of this expression is one, and M a/ 3 7 = otherwise. 

Maximum quantum value, Q m . We get a lower 
bound for the quantum value with measurement oper- 
ators given in Eq. ([7]) applied to components of a 3-qubit 
GHZ state, with Oj = 6® = 0% = and tp£ = <p% = 
ip^ — irfi/m. Using Eq. ([6]) it is clear that each nonzero 
Bell coefficient will be multiplied by the same value as 
itself, therefore, the quantum value will be equal to the 
sum of the absolute values of the Bell coefficients, that is 
with the no signalling limit, which is an upper bound for 
the quantum value. Hence it has the property of pseudo- 
telepathy [l5j]. From this it follows, that the quantum 
value will be nothing else than the number of nonzero 
Bell coefficients, which is actually m 2 . To see this, it is 
enough to note that however we slice up the m x m x m 
arrangement, each resulting 777 x m matrices will have 
exactly one nonzero number (plus or minus one) in each 
of its rows and columns. We can get such a row or col- 
umn by fixing two of the indices of M Q( g 7 . The sum of 
the indices we get this way are m consecutive numbers, 
exactly one of them will be divisible by m. Such a matrix 



will have m nonzero elements, the 777 slices together will 
contain m 2 such elements, therefore the quantum value 
and the no signalling limit will be Q m — m 2 . 

We now place an upj>er bound on the maximum of 
the biseparable value (B l m ), and then we prove that this 
bound can be saturated, that is, [B m = B m ). 

Upper bound on the biseparable quantum value, 

B^. A further property of the slices of the present m x 
777 x m arrangement is that they are modified circulant 
matrices like in the case of the previous family, which can 
be shown exactly the same way as we have shown there. 
To get the matrices relevant to the biseparable value, 
we have to add up the slices with different signs. If it is 
Alice who is not allowed to share entangled state with the 
others, this sum is = J2a=o AaMap^. Due to the 

property of the arrangement, for each matrix element, all 

terms of the sum but one will be zero. Therefore, each 
{A} 

entry of Ml will either be one or minus one. Moreover, 
this matrix will also be a modified circulant one, and its 
first line, which determines all the others, may contain 
any combination of plus and minus one values, depending 
on A a . Let a m _i_ 7 = J2™=o A a M a oy, that is the first 
line of written in opposite order. Then an upper 

bound for the biseparable value may be written as: 



B u 



E 

7=0 



cue 



(24) 



Let us consider the case of j — 0. Then what we get 
is the same as Eq. (fl~8|) but with A' = (which is ir- 
relevant), and a prefactor of m instead of m 2 /2. Then, 
according to Eq. (fT9)) , the result is mj sin(7r/2777). We 
will show that this actually is the upper bound, when- 
ever m is a power of two. As we have discussed earlier, 
exp (7717/777), which corresponds to j = 0, will point to- 
wards consecutive corners of a regular polygon of 2m 
edges on the complex plane while 7 takes all values be- 
tween zero and 777 — 1. If m is a power of two, then for 
any j, exp(iir(2j + 1)7/777) will point towards different 
corners for the different 7 values, moreover if one of them 
will point towards one corner, there will be none pointing 
towards the opposite corner. The reason is that (2j + 1)7 
is never divisible with m in this case. Choosing et 7 ap- 
propriately one can achieve that the terms to be added 
point towards consecutive corners, if taken in some order, 
which maximizes the absolute value of the sum. This is 
not true if 777 is divisible with an odd number. When 
2j + 1 is equal to this number, for 7 = m/(2j + 1) the 
value of exp(7T7(2j + 1)7/777) = — 1, which lies opposite 
to +1, the value for 7 = 0. In this case not all corners 
can be reached with appropriate choices of a 7 , and the 
other corners can be reached more than once, and B" n 
may be larger than what we have calculated. If m is 
odd, for 2j + 1 = 777 with a 7 = — l 7 we even reach the no 
signalling limit. 

Now we show that we can actually reach the value of 
B? n — 777/ sin(7r/2?77). 
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Lower bound on the biseparable quantum value, 

B\ n . The coefficients of the reduced Bell inequality are 
the elements of the modified circulant matrix M^ 7 whose 
entries in the first line are all +1. The appropriate Eu- 
clidean vectors C 7 are the same as the ones already de- 
fined in Eq. (|2"uT) . and analogously to Eq. (j2"Tj) we may 
write 



3=0 



7=0 



(25) 



For (3 = 0, Mo-y = 1, and we have to sum m unit vec- 
tors pointing towards consecutive corners of a polygon of 
2m sides, the usual formation, the length of the resulting 
vector is 1/ sin(7r/2m). For (3 = 1 only the last element 
of the row will be — 1. But if we change the sign of just 
the last vector of the formation, we get the same forma- 
tion rotated by an angle of w/2m. This formation will 
give the same result. The next line will give a formation 
rotated further by it /2m, and so on, therefore the result 
is raj sin(7r/2m). This is a lower bound for the bisepa- 
rable value, which is equal to the upper bound if m is a 
power of 2. In this case the ratio of the quantum and the 
biseparable limits is Q m /B m = msin(7r/2m), resulting in 
the threshold visibility V = B m /Q m = l/(msin(7r/2m)). 
This is the same visibility obtained under Eq. (|2"2l by 
means of the Bell inequality of section [TT1 

However, we would like to mention that this family of 
inequalities is more economical than our previous one. 
Namely, the number of joint measurements involved in 
Bell inequality ([5]) scales as m 3 , whereas the present Bell 
inequality defined by (|23|) consists of only m 2 joint mea- 
surements. Even for smaller number of measurements, 
the case which is more relevant to experiments, the dif- 
ference is not negligible: Inequality §5$ (or equivalently 
the Bancal et al. inequality [11|) gives the threshold vis- 
ibility V = 0.666, requiring 18 joint correlation terms. 
On the other hand, the inequality defined by (|2"5|) yields 



the lower threshold V = 0.653, using only 16 joint terms. 



IV. CONCLUSION 

In this paper we extended the three-party three-setting 
inequality of Bancal et al. [ll[ , which serves as a device- 
independent genuine tripartite entanglement witness, to 
an arbitrary number of settings. Our Bell inequalities 
(see Eq. ([5]) and Eq. (|2"3"|) for their definitions) can detect 
genuine tripartite entanglement in the noisy 3-qubit GHZ 
state with a visibility threshold of V = l/(msin(7r/2m)), 
where m denotes the number of settings per party. For 
m = 2, 3 our result recovers the threshold values corre- 
sponding to the Mermin inequality [l2| and the Bancal 
et al. inequality respectively. Numerical optimiza- 
tion suggests that these threshold values are optimal for 
m = 2 and m = 3. However, it is still an open question 
whether our family of inequalities (defined by Eq. ([5])) 
is optimal for any value of m. The optimality of these 
inequalities for m > 3 is supported by the fact that for m 
going to infinity the visibility V approaches the value of 
2/ir, achievable by a biseparable model simulating three- 
party correlators. Also, it would be desirable to general- 
ize our families either to more parties (here the method of 
Appendix C in Ref. [11] might be instructive) or to more 
outcomes. One may also wonder whether the inequalities 
presented in this work are optimal for important states 
different from the 3-qubit GHZ state. Furthermore, it 
would be of interest to find a Bell inequality, which is 
not the sum of three-party-correlators, giving a thresh- 
old visibility lower than 2/n for the noisy 3-qubit GHZ 
state. 
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